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ABSTRACT:  A  theory  i*  developed  for  the  effect  of  the 
earth's  rotation  on  its  spectrus.  Each  line  is  resolved 
by  rotation  into  a  multiplet  of  (2n-fl)  lines,  as  in  the 
Zeeman  effect,  A  theory  is  also  given  of  the  intensity 
distribution  in  the  lines  of  the  multiplet  for  the  case 
of  a  point-  source.  Good  agreement  is  obtained  between 
this  theory  and  the  doublets  observed  seismicaliy  and 
gravime tr i ca 1 ly . 
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Rotational  Multiplets  in  the  Spectrum 
of  the  Earth 

by 

C.  L.  Pekeris^  Z.  Alterman  and  H.  Jarosch 
Department  of  Applied  Mathematics,  The  Veizmann  Institute,  Rehovot,  Israel 


Abstract.  The  doublets  in  the  spectrum  of  the  free  oscillations  of  the  earth 
which  have  been  observed  on  the  gravimetric  (UCLA)  and  strain-meter  (Pasadena) 
records  of  the  great  Chilean  earthquake  of  May  22,  1960,  are  interpreted  as 
multiplets  arising  from  the  rotation  of  the  earth.  The  phenomenon  is  similar 
to  the  Zeeman  effect,  and  is  indeed  a  realization  of  the  mechanical  analogue 
from  which  Larmor  deduced  the  "Larmor  precession"  in  his  interpretation  of  the 
Zeeman  effect.  A  first-order  perturbation  calculation  yields  the  result  that 
the  degenerate  frequency  cro(n)  in  the  absence  of  rotation  is  resolved  by  a 
slow  rotation  into  (2n<t-l)  lines  an”*  given  by 

an"*  =  cro(n)  +  iirr(n)«  ,  -n  <  m  <  n  , 

where  u  denotes  the  angular  velocity  of  rotation  of  the  earth,  and  m  is 
the  azimuthal  number  of  the  wave  function.  r(n)  is  determinable  from  the 
zero-order  solution  in  the  case  of  spheroidal  oscillations,  and  is  equal  to 
[n(n-t-l)]~^  ,  in  the  case  of  purely  torsional  oscillations.  The  relative  inten¬ 
sities  within  tne  quintet  n  =  2  and  the  septet  n  =  3  have  been  determined 
for  an  observing  station  at  Los  Angeles,  on  the  assumption  of  an  explosive 
point-source  at  the  earthquake  focus  in  Chile.  The  strongest  lines  should  be 
the  pair  m  =  +1  for  n  =  2  ,  and  the  pair  m  =  +2  for  n  =  3  •  These  agree 
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in  s'^paration  with  the  pairs  observed  on  the  strain-meter  and  with  the  gravi¬ 
metric  pair  at  n  =  3  »  but  less  so  with  the  gravimetric  pair  at  n  =  2  , 

There  are  indications  in  the  strain-meter  spectrum  for  n  =  3  of  a  weaker 
line  at  m  =  0  ,  while  the  other  lines  are  theoretically  of  an  intensity  not 
exceeding  the  background  noise.  The  separation  in  the  observed  gravimetric 
doublet  for  the  first  overtone  of  n  =  3  agrees  with  the  interval  of  the 
strongest  pair  m  =  i2  .  The  intensities  of  the  lines  in  the  rotational 
multiplets  of  the  components  of  displacement  for  an  observing  station  at 
Palisades,  N.  Y.,  have  also  been  determined. 

Introduction. 

The  first  attempts  of  Zeeman  to  observe  an  effect  of  a  magnetic  field 
on  spectral  lines  led  to  negative  results.  He  resumed  the  experiment  in 
1894  when  he  read  in  Maxwell's^ sketch  of  Faraday's  life:  "Before  we  de¬ 
scribe  this  result  we  may  mention  that  in  1862  he  made  the  relation  between 
magnetism  and  light  the  subject  of  his  very  last  experimental  work.  He  en¬ 
deavored,  but  in  vain,  to  detect  any  change  in  the  lines  of  the  spectrum  of  a 
flame  when  the  flame  was  acted  on  by  a  powerful  magnet."  Zeeman  states^ 

"If  a  Faraday  thought  of  the  possibility  of  the  above-mentioned  relation,  per¬ 
haps  it  might  yet  be  worthwhile  to  try  the  experiment  again  with  the  excellent 
auxiliaries  of  spectroscopy  of  the  present  time,...."  The  discovery  of  the 
Zeeman  effect  that  followed  in  1896  came  at  a  time  when  the  bajic  concepts 
needed  for  its  interpretation  had  already  been  formulated  by  Lorentz^’^. 

Zeeman  goes  on  to  say^ that  "Professor  Lorentz,  to  whom  1  communicated  these 
considerations,  at  once  kindly  infcn.ied  me  of  the  manner  in  which,  according 
to  his  theory,  the  motion  of  an  ion  in  p  magnetic  field  is  to  be  calculated, 
and  pointed  out  to  me  that,  if  the  explanstion  following  from  his  theory  be 
true,  the  edges  of  the  lines  of  the  spectrum  ought  to  be  circularly  polarized." 
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The  splitting  of  the  original  frequency  Oq  by  a  magnetic  field  H  ,  which 
Zeeman  derives^ on  the  basis  of  the  Lorentz  theory,  is 


a 


00 


(1) 


for  the  case  where  the  lines  are  viewed  parallel  to  the  magnetic  field.  The 
two  lines  given  by  (1 )  should  be  circularly  polarized,  and  this  was  con¬ 
firmed  already  in  the  first  experiment  of  Zeeman^ 

The  relation  (1  )  had  been  derived  independently^*^  by  Larmor.  When 
news  of  Zeeman's  discovery  reached  him,  he  substituted  the  value  of  the  mass 
of  hydrogen  for  m  in  (1 )  and  concluded  that  the  effect  would  be  inappreci- 
'i.ble.  He  therefore  asked  Lodge  to  confirm  the  experiment.  This  Lodge  suc¬ 
ceeded  in  doing,  and  on  May  20,  1897>  he  demonstrated  the  effect  at  a  Royal 
Society  soiree.  Lodge  latei  published  a  disclaimer^ of  "any  intention  of 
trespassing  on  the  prerogative  of  the  discoverer*"  As  is  well  known,  Larmor 
proves  that  the  effect  of  a  uniform  magnetic  field  H  on  the  orbit  of  a 
particle  of  charge  e  and  mass  m  is  to  set  the  whole  orbit  into  a  pre- 
cession  around  the  direction  of  H  with  an  angular  velo''ity 


(2) 


provided  the  centrifugal  force  can  be  neglected  in  comparison  with  the 
Coriolis  force.  (ii|^  is  called  the  Larmor  precession  frequency.  Larmor 's 
theorem  follows  directly  from  equating  to  zero  the  sum  of  the  Coriolis  force 
2m  V  X  u  and  the  Lorentz  force  -(e/c)  v  x  H  .  Larmor  tken  shows,  as  Zee- 
man  had  done,  that  the  frequency  of  a  harmonic  oscillator  Oq  is  split  by 
a  rotation  u  into  the  two  lines  given  by  (1),  provided  (w/co)  «  1  . 
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While  the  Zeeman  effect  furnishes  ample  experimental  verification  of 
the  electromagnetic  part  of  Larmor's  analogy,  an  experimental  demonstra¬ 
tion  of  the  mechanical  counterpart  has  been  wanting.  A  striking  case  of 
the  splitting  of  the  natural  frequency  of  a  purely  mechanical  system  by  ro¬ 
tation  presented  itself  recently  when  the  records  of  the  great  Chilean  earth¬ 
quake  of  May  22,  1S60  were  analyzed.  This  earthquake  excited  the  natural  os¬ 
cillations  of  the  earth.  The  gravest  modes  n  =  2  (53»7  min)  and  n  =  3  (35*5  min) 
appeared  as  doublets  in  the  spectra  of  both  the  gravimetric  records^ *^and  of 
the  strain-meter  records^ 

The  free  oscillations  of  the  earth  are  governed  by  gravitational  and 
elastic  forces.  Taking  the  distribution  of  density  p(r)  and  of  the  elas¬ 
tic  constants  A(r)  and  p(r)  ,  as  inferred  primarily  from  seismic  as  well 
as  from  other  geophysical  data,  we  have  determined  the  spectrum  for  several 
proposed  models  of  the  earth^®^  / The  investigation  was  prompted 
originally  by  an  observation  made  by  Benioff^ ' of  a  57-minute  oscillation  on 
the  record  of  the  Kamchatka  earthquake  of  1932.  By  the  time  the  theory  was 
worked  out  in  detail  and  the  spectrum  thoroughly  investigated^  several  years 
had  elapsed;  and  a  cpiestion  began  to  emerge  with  annoying  persistence  as  to 
w^y  no  further  natural  oscillations  had  been  recorded  since  1932»  in  spite  of 
the  continuous  improvement  in  recording  facilities.  The  answer  to  this  ques¬ 
tion  came  at  the  meeting  of  the  International  Union  of  Geodesy  and  Geophysics 
held  at  Helsinki  in  July  I960,  when  Benioff-Press-and  Smith,  and  Ness-Harrison- 
and  Slichter  announced  that  they  had  identified  the  free  oscillations  of  the 
earth  in  the  strain  seismograms  (52  lines)  and  the  gravimetric  records  (49  lines), 
respectively,  of  the  great  Chilean  earthquake  of  May  22,  i960.  The  periods  de- 
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duced  from  spectral  analysis  of  the  strain-meter  records  agreed  with  the  gravi¬ 
metric  values  to  within  1  Vo  1  and  this  was  also  the  measure  of  agreement  with 
the  theoretical  spectrum^  for  the  Gutenberg  earth  model  and,  to  a  lesser  ex¬ 
tent,  with  the  spectrum  for  the  Sullen  B  model.  Free  modes  ranging  from  spheri¬ 
cal  harmonic  order  n=2  up  to  n  =  38  were  identified  seismically,  and  up  to 
ns^fl  gravimetrically. 

The  free  oscillations  fall  into  two  classes:  spheroidal,  with  non¬ 
vanishing  dilatation,  and  torsional.  In  the  latter  there  is  neither  dilatation 
nor  vertical  displacement,  so  that  they  produce  no  gravity  perturbation.  Actu¬ 
ally,  the  torsional  oscillations  were  identified  only  on  the  seismic  records. 

Press  and  Slichter  reported  that  the  n  =  2  and  -  3  lines  appear  as 
doublets  in  the  strainnneter  and  gravimetric  spectra.  It  was  then  suggested^ 
that  the  splitting  is  die  to  the  earth's  rotation.  This  conjecture  was  made  on 
the  basis  of  a  recollection  of  a  result  in  Lamb's  treatment^ of  the  effect 
of  rotation  on  the  free  gravitational  oscillation  of  a  circular  basin,  namely 
that  the  wave  advancing  in  the  direction  of  rotation  has  a  longer  period  than 
the  wave  going  in  the  opposite  direction.  A  perturbation  calculation,^"^ 
based  on  Lamb's  analysis,  showed  that  in  the  case  of  the  rotating  circular  ba¬ 
sin  the  frequency- interval  -  Oi)  in  the  doublet  should  be  of  the  order 

of  the  angular  velocity  of  rotation  u  .  Observationally,  the  quantity 
(oa  -ai)/u  ranged  from  0.7  to  1.1,  thus  lending  support  to  the  hypothesis  of 
t^'e  rotational  origin  of  the  observed  doubling  of  the  periods  of  free  oscilla¬ 


tions  of  the  earth. 
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In  the  following  section^  we  extend  our  previous  analysis  of  the  free  os¬ 
cillations  of  a  non-rotating  self-gravitating  elastic  earth  by  carrying  out  a 
first-order  perturbation  calculation  of  the  effect  of  a  slow  rotation  on  the 
frequency.  Denoting  by  subscript  zero  the  solution  for  the  case  of  no  rota¬ 
tion,  the  components  of  displacement  Uo,Vo,wo  in  a  spherical  system  of  co¬ 
ordinates  (r,d,<f))  and  the  perturbation  in  gravity  ,  are  given  by 


Uo  =  Uo(r)YnB(5»^)  » 


31 


Wo  = 


.  Vairiav, 


sind  d<f) 


Vo  =  V„(r>  gj 
i/iq  =  Po(r)Ynm 


(3) 


where 


Ynm(fl»<#>)  =  Pn^Ccose)  e^“^  , 


(4) 


iat 


and  a  factor  e*"'”  has  been  omitted.  The  frequencies  ao(n)  are  degenerate, 
and  do  not  depend  on  the  azimuthal  number  m  .  It  is  shown  below  that  the  in¬ 
troduction  of  a  slow  angular  rotation  u  («  oo)  removes  the  degeneracy,  each 
line  ao(n)  being  split  into  a  multiplet  of  (2n4-l)  lines  an”  given  by 

an"  =  ao(n)  +  mr(n)w  ,  -n  <  m  ^  n  (5) 

where 


pr* 

T(n)  =0 


(2UoVo  +  Vo*)dr 


pr*[Uo*  +  n(n+l)Vo*]dr 


(6) 


in  the  case  of  spheroidal  oscillations,  and 


'  (V) 

in  the  case  of  torsional  oscillations. 

Taking  the  earth  model  Bullen  B,  for  which  we  have  evaluated  the  func¬ 
tions  Uo(r)  and  Vo(r)  ,  we  find  from  (6)  that  t(2)  =  .395  and 
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t(3)  =  .183.  Using  these  in  (5)»  and  the  value  1,2~12  x  10’*  sec"^  of  the 
angular  velocity  of  rotation  of  the  earth  (o  ,  we  get  the  periods  Tn"  (=27r/an") 
of  the  quintet  for  n  =  2  and  of  the  septet  for  n  =  3  shown  in  Table  1 • 

The  question  arises  as  to  why  only  two  lines  out  of  a  possible  five  were  ob¬ 
served  in  the  case  n  -  2  ,  and  again  why  only  two  out  of  a  possible  seven 
lines  in  the  case  n  =  3  7 


Table  1_ 

The  periods  Tn®  for  model  Bullen  B, 


55.33 

54.50 

53.70 

52.92 

52.15 


35.99 

35.83 

35.67 

35.50 

35.34 

35.18 

35.03 


Observed 

Gravimetric 

Seismic 

min 

min 

54.98 

54.7 

52.80 

53.1 

35.87 

35.9 

35.24 

35.2 
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This  leads  us  to  an  investigation  of  the  relative  amplitudes  of  the  lines 
within  a  multiplet.  The  relative  amplitudes  depend  on  the  nature  of  the  source, 
its  geographical  location,  and  the  location  of  the  observing  station;  as  well 
as  on  the  nature  of  the  quantity  that  is  observed,  —  whether  a  component  of 
displacement,  or  of  strain,  or  the  perturbation  in  gravity.  We  assume  an 
explosive  compressional  point-source  at  the  earthquake  focus  in  Chile 
($0  s  128**),  and  make  use  of  the  results  of  an  earlier  investigation^*^  of  the 
relative  amplitudes  of  the  modes  n  s  2  and  n  «  3  »  which  such  a  source 
excites.  Were  the  source  located  on  either  of  the  poles,  there  would  be,  ac¬ 
cording  to  our  theory,  no  rotational  splitting,  because  there  would  be  no  longi¬ 
tudinal  (East-West)  component  of  motion  w  ,  and  m  in  (4)  and  (3)  would  be  zero. 
Although  on  the  assumption  of  a  point-source  there  is  no  longitudinal  motion 
around  the  axis  passing  through  Chile,  there  longitudinal  motion  around  the 
axis  of  rotation  of  the  earth,  and  this  gives  rise  to  rotational  splitting. 

In  the  case  of  the  gravimetric  measurements,  we  use  the  expressions  for 
the  amplitudes  within  the  multiplet  of  the  perturbation  of  gravity,  while  in 
the  case  of  the  strain-meter  we  use  the  appropriate  theoretical  amplitudes  of 
the  longitudinal  strain  in  a  horizontally-placed  long  bar. 

These  theoretical  results  are  compared  with  observations  in  Table  3  and 
in  Figures  1  and  2.  It  is  seen  that,  in  general,  the  doublets  observed  are 
those  lines  of  the  multiplet  which  are  strongest  theoretically,  and  that  the 
missing  lines  would  not  be  expected  to  stand  out  above  the  observed  noise  level. 
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2.  Theory  of  rotational  aultiplets  in  the  spectrua  of  the  earth. 

In  the  absence  of  rotatloni  the  analysis  of  the  free  oscillations  of  the 
earth  proceeds  by  first  determining  the  equilibrium  static  solution*  and  then 
superimposing  on  it  a  perturbation  velocity-field  v  which  is  controlled  by 


the  elastic  and  gravitational  restoring  forces. 

The  static  solution*  which  we 

shall  designate  by  the  subscript  xero*  is  spherically  symmetrical  and  is  gov¬ 
erned  by  the  equations 

(8) 

(9) 

tV  =  -  4sCpo  1 

(10) 

it  denoting  the  gravitational  potential  and  G  the  gravitational  constant. 

If  the  earth  now  rotates  with  angular  velocity  w  *  equation  (8)  still  retains 
the  form<<*^ 

^Po  »  PoVt  »  (1 1 ) 

where  t  is  the  seopotential .  comprising  the  gravitational  potential  $  de¬ 
fined  in  (10)  and  the  centrifugal  term* 

♦  »  p  +  ^  w*r*sin*®  *  (12) 

in  a  spherical  system  of  coordinates  (r*9*^}  •  The  solution  of  (11)  is 

Po  ■  Po(*)  > 

It  can  be  shown^*’^  that  the  surfaces  t  are  ellipsoids  of  revolution  whose 
maximum  ellipticity  is  at  the  surface  £  ■  a  *  where  to  a  close  approximation 

r  ■  o[i  -I-  (i/297)(ir  -  cos*P)]  .  (14) 
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The  geopotential  surfaces  are  therefore  nearly  spherical  throughout  the  volume 
of  the  earth,  to  within  1  part  in  300.  In  the  following,  we  shall  therefore 
neglect  the  ellipticity  of  the  geopotential  surfaces  and  shall  assume  that  po  $ 
Po  and  go  are  functions  of  £  only.  One  can  also  demonstrate  the  smallness 
of  the  effect  of  the  centrifugal  force  from  the  relation 

7**  =  +  2m*  »  -  4eGp  +  2m*  (15) 

which,  with  p  w  5  >  gives  (2m*/i|.«Gp)  «  (l/ifOO)  • 

The  rotation  does  affect  the  motion  through  the  "deflecting  force  of  the 
earth's  rotation",  or  the  so-called  Coriolis^ force  C  which,  per  unit 
volume,  is  given  by 

^  s  2poV  X  m  •  (l6) 

Let  u,£>v  denote  the  components  of  displacement  in  a  spherical  system 
of  coordinates  £,d,f^  •  The  equations  of  motion  and  of  the  perturbation  in 
the  gravity  field  are^*^*^’^ 


P  |t»  ■  •in®  If  ■  *  . 

(17) 

p  ||J  -  2up  cos9  If  “  S  * 

(18) 

p  1^  ♦  2mp  cos8  ll^^  2mp  sin8  |f  ■  T  * 

(19) 

♦  =  To  +  #  , 

7•^^  ■  iieG(pA  +  u^)  • 

(20) 

Here  p  denotes  the  unperturbed  density  po  >  A  is  the  divergence  of  the 
displacement,  and  a  dot  denotes  differentiation  with  respect  to  £  •  jt,  8 
and  T  are  linear  functions  of  the  displacements  and  of  the  components  of  the 
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strain  tensor  t  as  shown  in  the  Appendix*  The  time  does  not  appear  ex¬ 

plicitly  in  R,  S  and  T  * 

Having  solved^*^*^ for  the  functions  Uo(r)  and  Vo(r)  in  (3) 
and  for  the  corresponding  9o(n)  fo**  case  of  m  «  0  ,  we  now  proceed  to 
carry  out  a  first-order  perturbation  calculation  in  the  small  parameter  a  de¬ 
fined  by 

a  =  (w/oo)  «  1  •  (21 ) 

Let 

a  as  Oo  *  aoi  «  u  ■  uo  ♦  aui  ,  S  ■  So  4-  aSf  ,  etc*,  (22) 
then  the  sero-order  equations  are 

poo*Uo  +  Ro  =  0  ,  poro*Vo  +  So  =  0  ,  (23) 

P9o*Wo  +  To  =  0  ,  -  4fK»(pAo  +  Uo3)  *  0  •  (24) 

The  terms  proportional  to  a  yield  the  equations  for  the  determination  of 


ui  )  Vi  f  wi  and  a  : 

p9o*Ui  +  R|  X  .  2pooffiUo  -  2ip9o*  'in^  Wo  ,  (23) 

pao*Vt  +  S<  =  -  2p9o9iVo  -  Ripoo*  cosS  Wo  ,  (26) 

-I-  Ti  X  -  2paoaifio  2ip9o‘  cosS  Vo  *  2ip9o'  sinS  Uo  ,  (2?) 

7Vi  -  LsfG(pAi  Uf3)  *  0  •  (28) 


Let  a  star  signify  the  complex  conjugate*  It  then  follows  from  equations  (23)  - 
(27)  that 

(ufRo  -  utRi  +  vfSo  -  vf St  ♦  wfTo  -  wfTi )  x 

=  [2p9o9i(uou8  *  VoV$  Wow8)  *  2ipoo'sinS(utwo  -  Uewt)  + 

+  2ipao'cosd(vtwo  -  Vow8)]  *  (29) 
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It  is  shown  in  the  Appendix  that,  by  virtue  of  equations  (25)  -  (28)  and  the 

boundary  conditions, 
a  7r  2Tr 

Jr*dr^sin0d0^d^(utRo  -  uJRi  +  vfSo  -  vJSi  +  wfTo  -  wjTt )  =0  «  (30) 

0  0  0 

Hence,  by  (29),  the  volume-integral  over  the  sphere  of  the  r.h.s.  of  (29)  is 
zero.  This  furnishes  a  relation  for  detemining  in  terms  of  oo  and  the 

zero-order  solutions  defined  in  (3).  The  result,  derived  in  the  Appendix,  is 

Oi  =  moo  T(n)  ,  (31 ) 

with  r(n)  given  in  equations  (6)  and  (7)* 

It  follows  from  equations  (31 )  or  (5)  that  each  spectral  line  in  (3)  of 
non-vanishing  azimuthal  number  m  is  shifted  either  positively  or  negatively, 
depending  on  the  sign  of  m  .  It  is  onij  when  the  mode  excited  is  symmetrical 
about  the  axis  of  rotation  of  the  earth  that  the  rotation  has  no  effect. 

It  may  be  noted  that  the  free  oscillations  of  a  self-gravitating  liquid 
Maclaurin  ellipsoid,  for  which  the  multiplet-separation  has  been  determined^ 
for  any  finite  u  ,  approach  the  relation  (3)  in  the  limit  of  vanishing  u  , 
with  r(n)  =  l/n  .  The  separation  in  the  normal  Zeeman  effect  is  also  given^^*^ 
by  a  relation  (3),  with  r  =  1  ,  m  =  m^  ,  the  nmignetic  quantum  number,  and 
(u  s  ,  the  Larmor  precession  frequency. 

2.  Line  Intensities  in  Rotational  Multiplets. 

We  shall  now  determine  the  line  intensities  in  the  rotational  multiplets  of 
the  terrestrial  spectrum.  In  the  first  instance,  we  shall  treat  an  explosive 
compressional  point-source  at  the  earthquake  focus  in  Chile.  The  theory  can  be 
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generalized  to  more  complicated  types  of  source.  For  each  normal  mode>  the 
solutioh  of  equations  (23)  and  (24)  yields"^he  relative  values  of  Uo(a)  and 
Vo (a)  entering  in  (3)*  as  well  as  of  the  factor  P(a)  of  Ynm  in  the  per¬ 
turbation  of  gravity  at  the  surface  £  =  a  •  The  relative  excitation  of  the 
various  modes  by  a  compressional  point-source  has  been  determined^ for 
model  Bullen  B,  and  the  relevant  amplitudes  at  the  surface  are  shown  in 
Table  2.  These  were  determined  for  a  source  situated  at  the  surfacei  but  the 
relative  values  do  not  change  by  more  than  about  1  %  as  the  depth  of  focus  in¬ 
creases  to  200  km. 


Table  2 


Relative  excitation  of  modes  n  =  2  and  n  =  3  of 
model  Bullen  B  for  a  compressional  point-source  at  the  surface. 


n 

Uo(a) 

Vo  (a) 

^o(a) 

2 

1.00 

.0260 

.359 

3 

1.520 

00 

. 

1 

2.28 

Let  us  choose  a  spherical  system  of  coordinates  with  the  polar  axis  pass¬ 
ing  through  the  earthquake  focus  in  Chile  (38'*S,  73*3°W)»  and  let  B'  denote 
the  polar  distance  from  Chile.  Then  the  components  of  displacnent  u*  ,  v*  , 
w*  and  the  perturbation  of  gravity  at  the  surface  P*  are  given  by 


=  A(B') 

9 

CM 

1 

sinB' 

^  »<»•) '  0  . 

P'  =  C(B')  , 

(33) 
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where 

A(0’)  =  P2(co80')  e^^*^  +  1*520  P3(cos0')  +  .*•  ,  (34) 
B(fl’)  =  .0260  P*(cos6’)  -  .183  PsCcosS')  +  ...  ,  (35) 
C(e')  =  .359  Pa(co8e‘)  +  2.28  P3(cose')  +  ...  .  (36) 


Here  Oz  and  aj  denote  the  frequencies  for  n  =  2  and  n  =  3  respectively. 

In  order  to  determine  the  rotational  splitting  of  the  natural  oscilla¬ 
tions,  we  must  transform  the  above  quantities  to  the  geographical  coordinate 
system  (r,fl,<^)  referred  to  the  axis  of  rotation  of  the  earth.  6  denoting 
colatitude  and  <fi  longitude.  First  we  transform  A(0')  ,  B(d')  and  C(d') 
into  ^idt<t>)  ,  ^(Q$4>)  and  C(d,^)  respectively,  by  using  the  relation 


m=n 

(cose*)  =  Pn"'(cosfi)P„"(co8eo)e^(^“)  ,  (37) 


m=-n 


where  Pn"'(cos0)  «  Pn"*"  (cosd)  are  the  associate  Legendre  polynomials.  Here 
do  and  <Po  denote  the  colatitude  and  longitude  of  the  earthquake  focus  in 
Chile,  and  8  and  (p  the  colatitude  and  longitude  of  the  observing  station. 
The  components  of  displacement  u,v,w  and  of  the  perturbation  in  gravity  P 
in  the  geographical  system  of  coordinates  are  then  given  by 

U  =  A(»,#)  ,  V  =  I5  ,  (38) 

*  =  I? 

with  V  and  w  now  denoting  displacements  in  the  directions  of  South  and  East 
respectively. 

Thus  we  get  for  the  vertical  component  of  displacement  u  : 
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u  =  h{d,4>)  =  ,069  Pticoad)  -  .243  + 

+  .078  P8*(co80)[e^^^^“^^®^  +  + 

•/ 

t  e‘®'*^.517  P,(co89)  t  .134  P,’(co.e)[e‘(*^»)  +  ,-i(*-*o)]  . 

-  .073  Pj*(co.9)[e*(^^*»^  +  ,-i(2^-2*o)2  ^ 

+  .019  P,"(co8e)[«*(^^*^‘’^  +  e-l(3*-3*o)3'j  ^  .  (y3) 

J 

The  rotational  separation  of  the  multiplet  steins  entirely  from  the  dependence  on 
^  of  A(ef^)  ,  i(d,4>)  and  C{6,4>)  . 

The  Fourier  analysis  of  the  multiplets  for  n  =  2  and  n  s  3  was  made  by 
Benioff,  Press  and  Smith^^^  from  the  strain  seismograms  obtained  at  Isabella, 
California.  For  the  interpretation  of  their  results,  we  must  derive  an  expres¬ 
sion  for  the  strain  e  in  a  long  quartz  tube  anchored  to  the  ground  in  a  hori¬ 
zontal  position,  with  its  axis  making  an  angle  p  with  the  East  direction. 

We  have^'*^ 

e  =  (co8*/9  e^  ♦  ainpcoap  4-  sin*/9  e^^)  ,  (41 ) 

e^j  denoting  the  components  of  the  strain  tensor.  Using  equation  (3),  we  get 
for  each  partial  strain  en"*  corresponding  to  the  spherical  harmonic 
Y  =  Y„,(B,^)  , 

en*  =  UnY  +  Vn  [co8*/S(coW  ^  Y)  +  sin/5co8^(U  -  cotWY)  + 

+  sin'^ljl  ]  .  (1*2) 

Tills  leads  to 

ae  =  [A  +  cosVCcotW  ||  +  |^)  +  ^  sin^cos^  (ff  “  + 

4.  ain*p  ]  ,  (43) 

where  a  denotes  the  radius  of.  the  earth,  A  s  A(9,^)  ,  and  B  =  B(B,^)  . 
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We  note  that  this  expression  for  the  strain  along  a  horizontal  axis  contains 
the  term  A(d,^)  stemming  from  the  vertical  component  of  displacement  u  • 

Since  the  ratio  of  A/B  ,  as  given  in  (34)  and  (33) »  is  about  40  for  n  =  2  , 
and  about  8  for  n  =  3  «  it  follows  that  the  first  term  in  (43)  predominates; 
so  that  in  this  case  the  response  of  the  strain  seismograph  is  not  sensitive  to 
the  orientation  of  its  axis.  Indeed,  due  to  the  predominance  of  the  A-term,  the 
response  of  the  strain  seismograph  in  the  various  lines  of  the  muliiplet  re¬ 
sembles  the  relative  distribution  of  amplitude  in  the  vertical  component  of  dis¬ 
placement  u  ,  and  in  the  gravity  perturbation  ^  . 

Using  the  expressions  (38),  (39)  And  (43) »  we  have  computed  the  periods 
and  amplitudes  in  the  rotational  multiplets  of  the  Chilean  earthquake  for  ob¬ 
serving  stations  situated  at  Los  Angeles  and  at  Palisades,  N.  Y.  For  Los  An¬ 
geles,  the  gravimetric  amplitudes  were  computed  from  P  in  (39) »  and  the  seismic 
( so) 

strain  amplitudsa  Xrom  equation  (43) •  These  results  are  shown  in  Table  3  and  in 
Figures  1  and  2.  The  values  for  the  components  of  displacement  >  v  >  and  w 
in  Table  4  for  the  Palisades  station  are  based  on  equations  (38)  and  (39)*  The 
periods  in  the  multiplet  are  based  on  equation  (5),  with  the  values  of  r(n) 
given  in  Table  3* 

Before  proceeding  to  a  comparison  of  theory  with  observations,  we  shall 
discuss  the  Implications  of  our  assumption  of  an  impulsive  compressional  point- 
source.  To  begin  with,  the  tl  torsional  content  of  the  source  excited  tor¬ 
sional  free  modes  which  have  been  observed  seismically^  ,  and  which  lie  in  a 
different  part  of  the  spectrum.  Secmidly,  we  have  to  consider  the  evidence  put 
forward  by  Benioff,  Press  and  Smith^^^  to  the  effect  that  the  source  was  a  pro¬ 
gressive  rupture  proceeding  from  the  epicenter  southward  for  about  1000  km.  with 
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a  velocity  of  about  3  Vacate,  As  far  as  the  purely  compressional  content  is  con¬ 
cerned,  such  a  disturbance  could  be  represented  by  a  point-source  traveling  over 
a  distance  of  about  10°  of  latitude  and  lasting  sobk  5  minutes,  with  most  ct  the 
energy  probably  released  in  the  initial  burst  near  the  epicenter.  The  finite  time- 
extension,  as  against  an  assumed  8 -function  pulse,  nmy  affect  the  relative  excita¬ 
tion  of  the  entire  nodes  n*2  and  nK3  *  but  not  the  distribution  of  intensity 
within  the  nultiplet  of  each  mode,  fht  spatial  spreading  of  the  source  over  10° 
of  latitude  will  change  the  coefficients  for  the  relative  excitation  of  the  compo¬ 
nents  of  a  nultiplet,  such  as  is  given  in  (40),  to  a  degree  by  which  the  functions 
Pr’Ccos  00 )  vary  as  Co  goes  from  128®  to  138®.  For  n«  2  and  11  =  3  this 
variation  is  small,  and  is  of  the  sane  order  of  magnitude  as  the  a -terms  in  (22), 
which  we  have  neglected.  Since  a  =  (m/a)  is  around  l/2J|.,  we  may  expect  errors 
of  that  order  of  magnitude  fron  the  latter  source  in  the  aBg>litudes  given  in 
Tables  3  md  4. 

There  remains  to  be  considered  the  effect  of  torques  (SV)  whose  axis  has 
horizontal  components.  An  SV  torque  will  excite  spheroidal  oscillations  which 
are  not  synmetrical  about  an  axis  passing  through  the  source.  These  will  have  a 
different  distribution  of  intensity  in  the  rotational  nultiplets  from  that  of  a 
conprebsional  point-source. 
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Table  3 

Theoretical  perio4e-  T  and  amplitudes  of 
nultiplets  in  th«  spectrum  of  the  earth  excited  by  a  point- 
source  at  .Chile  28^)^  and  Abpirved  /at  Los  Angeles. 


Ampl. 

tude 

Observed 

periods 

n 

m 

T 

*min 

Gravimetric 

Strain 

Gravimetric^*^ 

Strain<^^ 

2 

-2 

55.33 

.115 

.302 

54.98 

2 

-1 

54.50 

.242 

.615 

54.7 

2 

0 

53.70 

.0015 

.0023 

2 

1 

52.92 

.242 

.615 

52.80 

53.1 

2 

2 

52.15 

.115 

.302 

3 

-3 

35.99 

.396 

.401 

3 

-2 

35.83 

1.257 

1.46 

35.87 

35.9 

3 

-1 

35.67 

.282 

.377 

3 

0 

35.50 

.624 

.873 

3 

1 

35.34 

.282 

.377 

3 

2 

35.18 

1.257 

1.46 

35.24 

35.2 

3 

3 

35.03 

.396 

.401 

Table  ^ 

Theoretical  periods  T  and  amplitudes  of  the  displacements  u  »  v  and  w 
of  nultiplets  in  the  spectrum  of  the  earth  excited  by  a  poTnt-source 
at  Chile  (9o  »  128^),  and  observed  at  Palisades,  N.  Y.  (9  b  49^). 


n 

m 

^min 

iL 

Vertical 

V 

South 

w 

East 

2 

-2 

55.33 

.265 

.0120 

.0182 

2 

-1 

54.50 

.721 

.0053 

.0248 

2 

0 

53.70 

.020 

.0053 

0 

2 

1 

52.92 

.721 

.0053 

.0248 

2 

2 

52.15 

.265 

.0120 

.0182 

3 

-3 

35.99 

.200 

.063 

.096 

3 

-2 

35.83 

.814 

.058 

_  .260 

3 

-1 

35.67 

.350 

.144 

.056 

3 

0 

35.50 

.288 

.162 

0 

3 

1 

35.34 

.350 

.144 

.056 

3 

2 

35.18 

.814 

.058 

.260 

3 

3 

35.03 

.200 

.063 

.096 
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Compari»on  of  theoretical  with  observed  multiplett  in  the  spectrum  of  the  earth. 

The  theory  presented  ^bove  predicts  that  every  frequency  of  spheroidal  oscil¬ 
lations  of  the  earth  of  order  n  is  split  by  rotation  into  a  multiplet  of  (2n>»-l) 
lines  whose  frequencies  are  given  by  equation  (5)»  m  denoting  the  asiauthal 
nunber  appearing  in  (3)  and  (k)»  Values  of  r(n)  are  shown  in  Table  5* 


Table  ^ 

Values  of  r(n)  in  equation  (6). 


_ fcl 

lien  B  mode 

(aO 

Gutenbersr  model  I 

i 

wng 

Ts 

wmtm 

T. 

■nm; 

I 

Overtone 

Ha 

min 

sec*' 

min 

2 

0 

1.931 

53.67 

.395 

1.957 

53.52 

.395 

2 

I 

4.235 

24.73 

.239 

4.306 

24.32 

.230 

2 

II 

6.764 

15.48 

.117 

6.912 

15.15 

.120 

3 

0 

2.953 

35.46 

.183 

2.964 

35.33 

.182 

3 

I 

5.8i*4 

17.92 

.213 

5.940 

17.63 

.207 

4 

0 

4.076 

25.69 

.099 

4.100 

25.54 

.098 

4 

I 

7.289 

14.37 

.198 

7.420 

14.11 

.192 

Using  these,  and  the  previously  determined  values^ of  the  central  frequencies 


(n)  ,  we  give  in  Table  3  the  periods  of  the  multiplet  lines  for  the  fundamental 
modes  of  spheroidal  oscillations  of  order  n«2  and  ns 3  i  In  the  case  of  model 
Bullen  B.  Also  shown  are  the  theoretical  amplitudes  of  the  lines  for  an  assumed 
compressional  point-source  at  Chile,  and  the  observed  periods* 

As  for  the  periods,  we  note  that  the  average  of  each  of  the  two  observed 
doublets  in  the  case  n«  2  is  53*9  nin,  as  against  the  theoretical  value  of  33*7 
for  model  Bullen  B,  In  the  case  ns3  >  the  average  period  of  the  observed  strain 
doublet  is  35*7  and  of  the  gravimetric  doublet  33*5^*  a*  compared  with  the  theore¬ 
tical  value  of  33*5  for  model  Bullen  B.  These  central  frequencies  depend  on  the 
model  assumed,  being  even  smaller  by  0.1 8  min  each  in  the  case  oC  the  Gutenberg 
model*  The  discrepancies  with  the  observed  values  disappear  for  the  Gutenberg 
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model  for  n>5  •  It  is  not  the  purpose  of  the  present  discussion  to  deal  with 
the  question  of  the  nodif  ications  which  would  be  required  in  the  structure  of  the 
model  in  order  to  eliminate  the  0*2  min  discrepancy  with  the  observed  central  peri¬ 
ods.  We  shall  rather  assume  that  the  whole  doublet  is  shifted  so  as  to  bring  the 
central  line  into  coincidence  with  the  observed  value,  end  we  shall  compare  the 
observed  intervals  between  the  lines  in  a  multiplet  with  the  theoretical  values 
for  the  intervals. 

A.  Cfaviifcrift.  Referring  to  Table  3*  in  the  case  n«2  the  central  line 
nsO  should  be  very  weak,  and  was  actually  not  observed.  The  theoretical  inter¬ 
val  between  the  strongest  jbb-i  and  msl  lines  is  1.6  min,  as  against  the  ob¬ 
served  interval  of  2.2  min.  Since,  however,  the  n«i2  lines  are  theoretically 
half  as  strong  as  the  ms±^  lines,  the  observed  lines  nay  represent  the  unresolved 
pairs  n  s  1  and  ns2  .  The  average  periods  of  the  pairs,  weighted  according  to 
their  amplitudes,  are  54*77  md  52.67  nin,  giving  an  interval  of  2.1  min,  whach  is 
close  to  the  observed  interval  of  2.2  nin. 

In  the  case  n«3  ,  the  strongest  lines  should  be  the  pair  £■  i2  ,  which  are 
very  close  to  the  observed  two  lines.  The  ■ «  0  line  should  be  half  as  stroi^,  and 
should  have  been  observed  if  the  noise  level  were  low  and  the  resolution  adequate* 

In  addition  to  the  doublets  in  the  fundoaental  n«  2  and  n»3  nodes,  ^e 
first  oveKone  of  the  ns  3  node  was  a  Iso  reported  as  a  gravitational  doubleti 
Ts17.88  and  17*68  min.  Theoretically,  the  relative  amplitudes  within  a  nnlti- 
plet  for  the  overtones  should  be  the  sane  as  for  the  fundanental,  as  given  in 
Table  3*  The  strongest  lines  in  all  the  overtones  of  iis  3  should  therefore  again 
be  the  pair  ns  i2  .  The  frequency  interval  is  Wri(3)  ■  6*20x  lO'^sec*^  ,  using 
the  value  of  .213  for  rx(3)  given  in  Table  5*  This  leads  to  a  period-interval 
for  the  pair  of  0.19  min,  which  agrees  with  the  observed  value  of  0.20  nin. 
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B.  Selwdc.  Two  doublets  were  rqixrted  by  Benioffi  Press  and  Snith^^^  trm  s 
Fourier  analysis  of  the  strain  seismogram  at  Isabella »  California*  In  Figure  1  is 
reproduced  their  spectral  intensity  curve  in  the  vicinity  of  the  pair  of 

54*7  and  33*1  min.  In  the  same  figure  are  shown  by  arrows  the  theoretical  separa¬ 
tions  and  amplitudes  of  the  5  lines  making  up  the  n«  2  quinteti  in  accordance 
with  the  results  given  in  Table  3*  The  central  £«  0  line  should  be  axtrssMly 
weak,  and  is  indeed  below  the  background  noise.  The  2  line  is  not  far  from 
the  observed  intensity  level*  On  the  other  hand,  the  a«-2  line  is  above  the 
observed  intensity  in  its  vicinity*  One  should  expect  a  rise  in  intensity  imme¬ 
diately  to  the  left  of  this  line,  which  is  outside  the  limit  of  the  figure* 

In  Figure  2  is  shown  the  observed  spectral  curve  for  n*  3  and  the  theore¬ 
tical  positions  and  amplitudes  of  the  £^*3  septet.  There  is  an  indication  of 
the  central  £bO  line,  and  its  theoretical  intensity  is  close  to  the  obeerved 
value.  The  pairs  are  close  to  the  genei^al  noise  level  in 

their  respective  neighborhoods* 

The  answer  to  the  question  posed  previously,  as  to  why  only  2  lines  out  of 
the  quintet  for  £■  2  have  been  observed  at  Los  Angeles,  is  that  the  observed 
doublets  are  the  theoretically  strongest  lines,  that  the  central  line  should  be 
week,  end  that  the  outermost  lines  would  not  be  expected  to  stand  out  above  the 
prevailing  noise  level.  In  the  case  of  the  septet  for  g"  3  ,  the  observed  doub¬ 
lets  are  again  the  theoretically  strongest}  there  is  actually  an  indicatien  of  a 
third  line  at  £«  0  in  the  strain  record,  while  the  missing  two  pairs 
and  £ai3  are  theoretically  again  Just  within  the  badiground  noise* 

It  would  be  of  interest  to  make  a  more  refined  Fourier  ena lysis  of  the  ob¬ 
served  spectral  curves  in  the  nultiplets,  and  also  to  determino  the  theoretical 
nultiplet  intensities  for  an  SV  torque-source*  Further  test  of  the  theory  could 
be  made  by  observations  at  other  latitudes,  liiere  the  intensity  distribution  within 
the  nultiplet  would  be  different.  We  are  at  the  threshold  of  a  new  sciencei 
terrestrial  spectroscopy. 
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APPENDIX  A 

1*  We  shall  first  prove  that,  by  virtue  of  equations  (23)  -  (28)  and  the 

boundary  conditions, 
a  tr  2r 

I  s  Jr*dr  JsinMd  Jd^(u*Ro  -  UoRi  +  v*So  -  vJSi  +  w*To  -  WoTi)  =  0  ,  (A1) 

0  0  0 

In  the  absence  of  rotation,  the  components  of  displacement  Uo  ,  Vo 
and  Wo  ,  and  the  perturbation  in  the  gravitational  potential  {fro  nre  given 
by 

Uo  =  Uo(r)Ynn(0,^)  ,  Vo  =  Vo(r)|y  Ynm  i  (A2) 

*0  =  Yn*  ,  {fro  =  Po(r)Ynm  ,  (A3) 

with 

Yn«(«,^)  «  P„*(cosD)e^^  .  (A4) 


We  have<*^*^*J 
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'#  =  rnsj  5?  *  F  ♦  F  '  'r*  =  -  F  *  F  fj 

1 

(A9) 

1  3u  8v  w 

*r^  ”  r  sinW  Jr  r  * 

(A10) 

'  F  n  -  F  ♦  F1IF?  H  • 

(A11) 

^  *  1?  +  ^  ♦  r  III.*  ♦  rUFF  5?  * 

(A12) 

We  now  substitute  froa  (A5)>  (a6)  and  (A?)  into  (A1)  and  carry  out 
the  indicated  integration.  Every  term  which  involves  derivatives  of 
Kf  Ht  L  and  of  the  components  of  the  strain  tensor  ,  we  integrate 
out  partially.  The  integrated  quantities  vanish  because  of 

1 )  the  boundary  conditions  at  the  surface  r  ~  t 

AA  +  2M(8H/ar)  •*  Cp#  *  «  0  ,  (AI3) 

2)  the  regularity  of  the  and  of  A  at  the  poles>  and 

3)  the  periodicity  condition  in  ^  •  It  is  found  that  the  terms  in  the 
integrand  multiplying  K  ,  n  ,  and  gp  each  cancel  out,  leading  to 

I  .  «?  U’*  H’-  v{  U’)  + 

where  jdr  denotes  the  volume-integral  in  (A1). 

For  the  proof  of  the  vanishing  of  (AI4),  we  make  use  of  the  equa¬ 
tions  for  the  gravity  potential} 


yVo  -  4»c(pAj  +  ujp)  *  0  ,  (A15) 

-  4sC(pAf  +  uf^)  «  0  .  (A16) 

Multiplying  (AI5)  by  and  (AI6)  by  -^0  $  adding  and  integrating 

over  the  volume  of  the  sphere,  we  get 

w  2w  ^  • 

Jsinddd  jd^*j^i^|jf®-  i^puS^  -  4sCpuT^  ■ 

•  If’-  «!  IS:'*  ?(’'?  W-  -i  H')  *  Hr-  -S  IS’|  •  ("7) 
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Now,  while  Uo  and  {fro  in  (A2)  and  (A3)  are  each  represented  by  a  single 


spherical  harmonic  of  order  n  ,  the 

tions  (23)  and  (28)  respectively  are 

solutions  for  Ui  and  {fri 

of  more  general  type 

from  equa- 

ui  ■  V**!*'^  » 

uU)  =  uU)(r)Yfc,(g,0)  , 

(AI8) 

Jt 

-  y#;*>  . 

^U)  „  pU>(r)Yfc,(g,^)  , 

(A15) 

K 

Uo  ■  uj^"^  , 

{fro  «  {fr‘,"^ 

(Aao) 

The  boundary  condition  at  the  surface  for  the  gravitational  potential 

isC  •),<•) 

-  4sGpu<‘'>  =  - ,  r  =  a  ,  (A21) 


for  each  k  •  Hence 

itOirnSj  -  i  )/,“>*•<•'  - 

-  .  (A22) 

The  surface-integral  of  (A22)  vanishes  because  of  the  orthogonality  of 
the  spherical  hamonics  when  n  /  k  ,  and  because  of  identical  vanishing 
tor  n  s  k  •  Hence,  the  volume-integral  in  (AI7)  vanishes,  and  with  it 
^  ,  as  given  by  (A14)« 

The  above  derivation  holds  for  both  spheroidal  and  torsional  oscil¬ 
lations,  since  all  the  boundary  conditions  are  satisfied  in  both  cases* 

In  the  case  of  the  torsional  oscillations,  some  of  the  boundary  condi¬ 
tions  are  satisfied  identically  because  for  them 

(A23) 


Uo  B  Ao  B  To  ■  {fro  0 
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2*  Evaluation  of  r(n)  . 


In  the  case  of  spheroidal  oscillations,  uo,  Vq  and  wo  are  given 
by  (A2}  and  (A3),  and  the  voluae-integral  of  the  r»h*s.  of  (29)  leads  to 


2or#ai  j drpj^Ujp*  +  \ 

•  (di)*  ♦  ]■  +  COM  p||  v; 

.(*») 

where 

P(®)  ■  Pj[(cosW)  * 

(A25) 

We  have 

*  . /.low,  p*  -  ^  , 

0 

(«<) 

w 

B  s  d$  =  n(n+1  )A  -  n*C  , 

0 

(A27) 

0 

(*28) 

> 

w 

nJcosB  ®  ^  A  * 

(A29) 

Using  these  in  (A24>),  we  are  led  to  ecpiations  (3)  and  (4)* 

In  the  case  of  torsional  oscillations,  we  have^*^*^*^ 

6 

o 

H 

o 

(A30) 

and  the  volume- integral  of  the  r*h*s*  of  (29)  leads  to 

2aoa,  jdrp  Vj 

(S)*  ♦  ii5*«  ■’*  "  • 

(A31) 

Fron  this,  relation  (?)  follows* 


MAGNITUDE  F(cii) 


m 


FREQUENCY  IN  CYCLES  PER  MINUTE 

FIGURE  I  . 

n  •  2  .  Spectral  intensity  of  Chilean  earthquake  observed 
at  Isabella  by  Benioff,  Press  and  Smith  .  Arrows  show 
theoretical  positions  and  amplitudes  of  multipict  for 
a  comprcssional  point  source  . 


amplitudes  of  mUltiplet  for  o  compressionol  point 
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